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ABSTRACT 


This  report  describes  22  separate  investigations  in  the  field  of  seismic 
discrimination.  These  are  grouped  as  follows:  Discrimination  and  evasion 
research  (5  contributions),  method  for  the  determination  of  focal  depth 
(4  contributions),  global  seismic  network  studies  (4  contributions),  studies 
of  earth  heterogeneity  (4  contributions),  source  mechanism  studies  (3  con- 
tributions), and  investigations  of  global  seismic  activity  (2  contributions). 
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SUMMARY 


This  is  the  twenty-st  annual  Technical  Summary  report  describing  the  activities 

of  Lincoln  Laborato^*  ne  field  of  seismic  r’iscr;rr  nation.  These  activities  involve 

research  into  the  fundamental  seismologicai  jm  associated  with  the  detection,  location, 

and  identificatijn  of  earthquikes  and  nuclear  explosions.  We  are  also  concerned  with  the  devel- 
opment of  methods  for  the  handling  and  analysis  of  large  quantities  of  global  seismic  data,  both 
from  the  point  of  view  of  data  management  system  design,  and  also  to  facilitate  the  optimum  ex- 
traction of  scientific  information  from  high-  quality  digital  data. 

Research  into  topics  directly  related  to  seismic  discrimination  and  evasion  has  concentrated 
on  two  areas.  First,  some  new  approaches  to  the  determination  of  surface  wave  magnitude  are 
described.  It  is  shown  that  useful  information  about  the  shape  of  the  group  velocity  curve  can 
be  deduced  directly  from  the  seismogram,  and  this  leads  to  a much  more  satisfactory  path  cor- 
rection than  those  commonly  employed.  Second,  we  are  exploring  multivariate  techniques  that 
will  enable  a much  more  complete  parametric  representation  of  the  seismogram  to  be  applied 
to  the  problems  of  discrimination  and  evasion.  Current  research  includes  the  simultaneous  ap- 
plication of  multiple  discriminants,  discrimination  between  different  source  mechanisms,  and 
an  approach  to  the  multiple  shot  evasion  problem. 

A number  of  studies  of  methods  for  the  determination  of  focal  depth  are  described.  Very 
promising  results  are  shown  for  the  identification  of  the  surface  reflection  pP  from  explosions 
using  maximum  entropy  cepstral  analysis.  This  technique  has  the  potential  for  becoming  a 
powerful  tool  in  any  seLsmic  monitoring  scheme.  Continued  study  of  the  effect  of  depth  on  sur- 
face wave  spectra  shows  a strong  influence  of  crustal  structure,  enough  to  make  the  method 
somewhat  suspect.  Application  of  measurements  of  S-P  times  is  discussed,  and,  although  the 
method  is  potentially  powerful,  there  is  clearly  considerable  difficulty  in  identifying  the  H ar- 
rival in  many  instances.  Another  study  attempts  to  model  the  characteristics  of  the  proposed 
DARPA  seismic  network.  It  is  shown  that  this  network  will  have  some  difficulty  in  resolving 
accurate  focal  depths  for  many  events  in  Eurasia. 

We  are  continuing  our  work  on  the  development  of  a management  structure  for  data  and  in- 
formation from  the  proposed  DARPA  seismic  network.  Application  of  the  NLS  system,  devel- 
oped by  Stanford  Research  Institute,  appears  pror  ising,  and  experiments  in  the  storage  id 
retrieval  of  information  via  the  ARPANET  are  being  carried  out.  Other  network  studies  include 
a continuation  of  our  eifort  during  the  International  Seismic  Month,  with  emphasis  on  the  char- 
acteristics of  a 32-station  network.  With  careful  data  analysis,  it  appears  that  such  a network 
is  almost  as  good  as  the  whole  set  of  ISM  stations  for  the  production  cf  a global  earthquake 
bulletin. 

Our  studies  of  the  nature  of  inhomogeneity  in  the  earth  are  becoming  heavily  involved  in 
scattering  processes,  particularly  in  the  mantle.  Several  investigations  described  in  this  report 
indicate  that  these  processes  may  play  a dominant  role  in  the  absorption  and  conversion  of  energy 
in  short-period  body  waves.  The  presence  of  anisotropy  in  the  scattering  medium  is  postulated, 
and  is  shown  to  constitute  a potentially  important  effect 

In  the  field  of  seismic  source  mechanisms,  we  describe  se'  ?ral  studies  aimed  at  the  inter- 
pretation of  source  parameters  from  observations  of  long-period  body  and  surface  waves.  One 
study  indicates  the  possibility  of  determining  gross  source  characteristics  from  long-period 
P-waves.  Another  describes  a method  for  the  estimation  of  source  parameters  from  the  am- 
plitude and  phase  spectra  of  surface  waves  from  pairs  of  earthquakes. 
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We  have  also  initiated  a study  of  the  time  vaiiation  of  seismic  activity,  hoping  to  clarify 
problems  rrlated  to  network  detection  and  hide -in-earthquake  evasion.  Evidence  is  shown  that 
the  level  of  g.obal  seismic  activity  shows  rapid  fluctuations  in  level  that  are  largely  independent 
of  magnitude.  Viewed  as  a time  series,  the  activity  variations  appear  to  contain  some  signifi- 
cant periodicities,  though  the  cause  of  these  has  yet  to  be  established. 

In  anticipation  that  digital  seismic  data  soon  will  be  available  in  quantity  on  the  ARPANET, 
we  are  continuing  to  develop  a facility  to  enable  us  to  access  and  manipulate  these  data.  We 
are  presently  in  the  process  of  developing  the  necessary  software,  and  we  anticipate  being  able 
to  access  this  nojv  information  source  when  it  becomes  available. 


M.  A.  Chinnery 
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SEISMIC  DISCRIMINATION 


I.  DISCRIMINATION  AND  EvASION 


A.  Mg:mb  FOR  EARTHQUAKES  WITH  COMPLEX  SHORT-PERIOD  WAVEFORMS 

It  has  teen  shown1  that,  by  firing  an  ordered  sequence  of  explo.  ms,  the  complexity  and 
Mg:mb  ratio  of  an  earthquake  may  be  mimicked.  It  is  reasonable  to  conjecture,  however,  that 
the  complexity  of  the  short-period  waveform  and  the  level  of  long-per  od  energy  contained  in 
the  surface  waves  are  both  due  to  the  source  mechanism.  The  rest’*  ,'ould  be  that  earthquakes 
of  a given  Mg:mb  ratio  would  have  a characteristic  short-period  waveform  as  well.  To  investi- 
gate this  hypothesis,  earthquakes  from  the  linear  region  of  shallow  focu.s  seismicity  just  north 
of  the  Hindu  Kush  wrere  divided  into  two  classes  based  on  the  emergent  or  impulsive  character- 
istic of  the  first  15  sec  of  their  short-period  teleseismic  P-waveforms,  and  their  Mg:m^  ratios 
were  computed.  Figure  1-1  shows  typical  waveforms  of  the  two  classes  and  a synthetic  seismo- 
gram produced  by  the  superposition  of  a P-waveform  of  an  explosion.  Dae  to  the  constructive 
interference  of  the  surface  waves  and  the  use  of  an  initially  small  explosion,  .he  Mg  will  increase 
by  approximately  0.5  and  '.he  will  decrease  by  approximately  0.5  for  the  multiple  explosion 
compared  with  that  of  a single  explosion.  Such  behavior  is  illustrated  in  Fig.  1-2  where  the  box 
labeled  A represents  the  spread  of  Mg:mb  at  - 5.7  5 for  presumed  Eastern  Kazakh  explosions^ 
and  the  box  labeled  B represents  the  expected  M :m.  region  in  the  earthquake  population  that 
the  multiple  explosion  would  fall  into.  The  points  marked  with  an  X arc  for  Class  2 (impulsive) 

earthquakes,  while  those  marked  with  an  O represent  Class  1 (emergent)  earthquakes.  The 
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trend  line  for  M„:mb  for  this  region  shown  in  the  figure  is  also  the  line  that  separates  the  two 
classes  of  short-period  waveforms.  While  the  multiple  event  Mg:mb  is  earthquake -like,  its 
waveform  (Class  1)  is  not  characteristic  of  events  with  comparable  Mg: m^  values. 

The  purpose  of  this  study  was  not  to  show  that  rnultiple-event  evasion  is  not  possible,  but 


rather  that  under  any  evasion  scheme  the  total  character  of  both  the  long-  and  short-period 


seismic  energy  must  be  considered. 


T.  E.  Landers 


R.  E.  Needham 


B.  SURFACE  WAVE  MAGNITUDE  '’ATI!  CORRECTION 

During  the  computation  of  surface  wave  magnitudes  (Mg)  for  the  International  Seismic 

Month  (ISM)  t.'vents,  it  became  clear  that  appropriate  path  and  period  correction^  were  necessary 

when  using  globally  distributed  stations  with  various  response  characteristics.  In  the  ISM  long- 

period  study,  M path  corrections  given  as  a function  of  period  by  Marshall  and  Basham^  were 
s 

applied.  Although  these  correclions  represent  a signiiicant  advance  in  reducing  the  variance 

in  an  M determ inaL  m from  a network  of  stations,  they  are  cumbersome  to  apply  in  any  auto- 
s 

matic  promising  • lemc.  The  difficulties  lie  in  determining  the  nature  of  the  path  of  any  given 
surface  wave  ray  and  the  magnitude  of  the  path  correction  along  paths  not  studied  by  Marshall 
and  Basham.  As  tnese  authors  correctly  point  out,  the  path  correction,  as  a function  of  period, 
depend,  chiefly  upon  the  slope  of  the  group  velocity  curve  between  the  epicenter  and  each  re  - 
cording  station.  The  suggestion  is  made  here  that,  in  some  cases,  this  slope  can  be  directly 
from  the  seismogram  and  a path  correction,  which  is  in  a sense  optimum,  may  be  computed  for 
each  individual  amplitude  measurement. 
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It  has  been  shown  by  Ben-Menahem  and  Jeffreys  and  others  that  the  amplitude  of  funda- 
mental mode  surface  waves  of  a given  frequency  (w  ) may  be  approximated  ir.  the  time  domain 

by 


A,'-A!3  /V,Uo)exp[,-%t/2Q)(TaT),/Z] 


(sin  A) 


(1-1) 


where  A is  the  epieentral  distance  in  degrees,  Lq  is  the  source  radiation  effect,  A(u,'o)  is  the 
source  spectral  amplitude  at  Q is  a dimensionless  attenuation  parameter,  and  is  the 
second  derivative  of  with  respect  to  wave  number  k at  u>q.  In  Eq.  (I-i),  the  parameters 
that  determine  the  source  strength  or  magnitude  at  a given  azimuth  and  frequency  are  and 
A(u>o).  The  other  terms  are  path  effects  for  which  we  would  like  to  correct.  If  we  assume  that 
the  epieentral  distance,  the  travel  time,  and  the  surface  wave  a.-enuaticn  parameter  are  known, 
then  the  only  remaining  correction,  the  path  correction,  is  proportional  to  (tl  d2u>/ak2  | o)*/2. 

It  may  be  shown  that  this  factor,  in  terms  of  quantities  that  can  be  directly  measured  on  the 
seismogram,  may  be  approximated  by 


V/2  * A ±°  (M.,1/2 

L2  . / t 6T 


Here,  Tq  is  the  period  (in  seconds)  of  the  particular  phase  where  magnitude  is  being  measured, 
t is  travel  time  (in  seconds),  5t  is  the  difference  in  arrival  time  of  the  preceding  and  succeeding 
phases,  and  6T  is  the  difference  in  the  period  of  these  phases.  Again,  the  latter  two  differences 
may,  in  some  cas*  ?,  be  measured  directly  on  the  seismogram  and  they  allow  the  reader  to 
account  for  the  path  dispersion  effects  in  a manner  that  is  relevant  to  that  seismogram. 

An  example  of  the  method  is  shown  in  Fig.  1-3.  The  data  trace  in  this  figure  represents 
the  vertical  component  of  the  Rayleigh  wave  recorded  at  the  World-Wide  Standard  Seismograph 
Network  (WWSSN)  station  at  Valent  ia.  Ireland  from  a Novaya  Zemlya  explosion.  The  triangles 
on  the  lower  plot  represent  Mg  units  measured  ly  the  conventional  term  log^0(A/T)  at  the  num- 
bered phases  of  the  seismogram.  The  open  circles  represent  the  term  log^A)  plus  a path 
correction  computed  as  suggested  above.  In  both  cases,  the  measurements  of  amplitude  (A) 
were  corrected  for  the  frequency  response  of  the  WWSSN  long-period  instruments.  It  may  be 
seen  in  Fig.  1-3  that,  in  this  case,  the  suggested  method  for  path  correction  gives  a more  stable 
Mg  measurement  over  a wide-period  range  than  the  conventional  method.  A test  of  the  suggested 
method  on  a large  data  population  is  being  undertaken.  The  n thod  will  not  apply  if  multipath 

propagation  is  evident  on  the  long-period  seismogram,  but  a few  M measurement,  made  by 

s 

the  new  method  it  ay  prove  more  reliable  *.han  many  measurements  made  on  complicated 
seismograms. 

An  alternative  method  for  estimating  the  slope  of  the  group  velocity  curve  is  described  in 
the  following  contribution.  . R .... 


C.  AUTOMATIC  AND  CONTINUOUS  MEASUREMENT 
OF  AMPLITUDE  AND  FREQUENCY 

As  reported  elsewhere  in  this  SATS,  a new  way  of  measuring  surface  wave  magnitudes  is 
being  developed.  This  method  reouires  that  the  slope  of  the  group  velocity  curve  as  well  as  the 
amplitude  and  frequency  of  the  arriving  signal  be  measured.  In  fact,  if  instantaneous  frequency 
and  amplitude  can  be  measured  as  a function  of  time  for  a surface  wave,  the  method  allows  one 
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to  measure  Mg  over  the  complete  range  of  arriving  frequencies.  We  have  started  to  investigate 
a method  by  which  such  measurements  can  be  made  automatically  by  a computer. 

A dispersed  Rayleign  wave  can  be  approximated  by 

s(t)  = A(t)  cos  [9(t)l 

where  A(t)  is  the  envelope  cf  the  signal  and 
w(t)  = ^ O(t) 

is  the  instantaneous  radian  frequency.  If  neither  w(t)  n ,r  A(t)  change  tco  quickly  with  time* 
there  is  a way  to  extract  A(t)  and  w(t)  which  is  quite  convenient  to  realize  automatically  on  a 
computer.  Given  w(t),  A(t)  is  the  amplitude  required  to  measure  Mg  at  that  frequency.  Since 
we  have  w(t),  the  origin  time,  and  the  distance  of  the  event,  we  can  also  obtain  the  required 
group  velocity  curve.  Thus,  the  time  functions  w(t)  and  A(t)  will  allow  us  to  obtain  a time  func- 
tion M (t).  The  value  of  'his  time  function  at  time  t will  be  an  M value  for  frequency  w(t)/2n. 
s s 

The  method  uses  the  complex  analytic  representation 
s(t)  = s(t)  + is(t) 

where  s(t)  is  the  Hilbert  transform  of  the  real  signal  s(t).  The  Hilbert  transform  can  be  realized 
:ther  by  a time  domain  filter  or  by  application  of  a 90°  phase  shift  in  the  frequency  domain. 

Both  methods  are  being  tried,  although  the  data  shown  below  were  obtained  using  the  frequency 
domain  approach.  The  functions  A(t)  and  w(t)  are  obtained  from  3(t)  by 

A(t>  = ^(t)  s*(t) 

and 

w(t)  = ^ Im  {In  s(t)} 

v/here  Im  denotes  taking  the  imaginary  part  of  the  complex  logarithm. 

Figure  1-4  shows  a seismogram,  the  envelope  A(t),  and  the  frequency  w{t)/z7r  obtained 
from  it.  Time  differentiation  to  obtain  w(t)  is  just  a scaled  first  difference  of  adjacent  data 
l .nts,  and  we  plan  to  consider  better  algorithms  which  will  not  give  such  a noisy  output.  The 
very  top  and  bottom  tra  s are  A (t)  and  w(t}/2.  filtared  by  a low- pass  filter  with  a corner  fre- 
quency of  0.02  Hz. 

This  effort  to  automatically  and  continuously  measure  amplitude  and  frequency  is  just  getting 
under  way.  The  initial  results  are  encouraging.  We  plan  now  to  improve  our  programs  and 
combine  them  with  algorithms  for  measuring  Mg.  The  computer  load  to  extract  A(t),  w(t)  is 
not  excessive  and  is  about  equivalent  to  filtering  with  a 30-point  convolutic nal  filter. 

R.  T.  Lacoss 

D.  APPLICATION  OF  PATTERN  RECOGNITION  TO  SEISMIC  DISCRIMINATION 

In  the  past  four  years,  there  has  been  a revival  of  interest  in  short-period  discrimination 

of  explosions  and  earthquakes,  particularly  at  low  magnitudes  for  which  surface  wave  measure- 
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ments  cannot  be  made.  Anglin  demonstrated  that  explosions  and  shallow  earthquakes  in 
Eurasie,  recorded  at  the  Yellowknife  array,  could  be  separated  on  a plot  of  complexity  vs  third 
moment  of  frequency,  whereas  if  either  parameter  were  used  alone,  there  was  considerable 
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overlap  betwe-  n the  explosion  and  earthquake  populations.  Further  studies^”7  verified  the  value 
of  multiple  discriminants  and  even  multisite  studies,  although  it  has  become  clear  that  careful 
regional  studies  of  spectral  parameters  must  be  made. 

We  are  attempting  to  apply  some  techniques  of  pattern  recognition  to  the  discrimination 
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proolem.  The  method  of  Fukunaga  and  Koontz  is  promising  because  it  is  designed  to  extract 
discriminating  features  for  two  classes  of  data.  Briefly,  the  technique  assumes  that  n param- 
eters have  been  measured  for  each  event  of  two  distinct  classes.  For  each  event,  the  various 
parameters  are  elements  x.  of  a data  vector 

X = col  (xrx2  . ..,xn)  . (1-2) 

For  each  class,  the  autocorrelation  matrix 


R.  = E.{XXT}  , i = 1,  2 


(1-3) 


is  estimated,  where  Ei  is  the  expectation  of  Class  i data.  A "mixture"  Rq  of  the  two  classes 
is  defined  as  tne  sum  of  the  two  autocorrelation  matrices,  i.e.. 


Ro  = R1  + R2 


(1-4) 


The  next  step  is  to  determine  a coordinate  transformation  F which  diagonalizes  and  whitens  R . 
Applying  this  to  Rq,  we  o!  tain 

I = S4  + S2  (1-5) 


whe  re 


S.  = PR.PT  , i = 1, 2 . (1-6) 

From  Eq.  (1-5),  it  can  be  shown  that  and  S,  have  the  same  set  of  eigenvectors,  and  for  each 
eigenvector  the  eigenvalues  of  Sj  and  are  related  by 

(1-7) 

where  1 >X^  > 0.  If  is  the  largest  eigenvalue,  e.g.,  nearly  unity,  then  i s almost  0,  so 

that  the  eigenvector  with  these  eigenvalues  is  the  most  important  feature  or  basis  vector  of 
Class  1 data,  while  simultaneously  being  the  least-significant  feature  of  Class  2 data. 

Projecting  the  data  vectors  of  each  class  onto  the  eigenvectors  with  eigenvalues  near  1 or 
0 should  give  the  best  separation  of  the  two  classes.  By  assuming  a normal  distribution  of 
each  data  class  as  projected  onto  these  etgemc  tors,  the  crossover  point  and  mis  identification 
probabilities  can  be  determined.  That  eigenvector  which  gives  the  smallest  m ^identification 
probability  is  therefore  the  direction  along  which  the  class  separation  is  best. 

To  illustrate  this  technique,  an  experiment  was  run  using  11  Central  Asian  earthquakes  and 
12  presumed  explosions  in  Eastern  Kazakh.  The  earthquakes,  Leken  to  be  the  first  class  of 
events,  are  from  the  Kirgiz-Sinkiang  border  region,  Tadzhik-Sinkiang  border  region,  End 
southern  Sinkiang,  all  within  about  10°  distance  from  the  concentration  of  Easterr  azakh  events, 
which  form  Class  2 data.  The  data  used  were  recorded  at  NORSAR,  and  the  parameters  in- 
vestigated were  m,  . lf>g  third  moment  of  the  spectrum  (LTM),  complexity  (CMPLX),  and  M . 
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The  significant  she  .-period  parameters  are  LTM  and  complexity,  which  are  plotted 

g 

in  Fig.  1-5  for  the  two  classes  of  data.  The  LTM  values  were  obtained  by  Noponen  by 
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spectraforming  the  NORSAR  subarray  beams  of  each  event.  This  method  of  estimating  short- 
period-event  spectra  is  designed  to  preserve  the  high  frequencies  usually  lost  in  conventional 
beamforming.  The  complexity  values  were  obtained  using  the  quadratic  complexity  formula  of 
Kelly ^ with  an  adcitional  correction  for  noise  power. 

The  separation  of  the  two  classes  is  very  clear,  although  either  discrimination  parameter 
alone  does  not  give  perfect  separation.  Figure  1-6  shows  an  M vs  m.  plot  for  the  events  of 
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Fig.  1-5  which  had  NORSAR  Mg  values  measured.  As  in  Fig.  1-5,  the  separation  is  perfect. 

Using  these  data,  the  Fukunaga-koontz  metnod  described  a novo  was  applied  to  the  two 
classes  of  events.  The  results  are  summarized  in  Table  I-i.  In  each  case,  the  parameters 
used  in  the  data  vectors  are  marked  by  x's. 


TABLE  1-1 

PATTERN  RECOGNITION  USING  CENTRAL  ASIAN  EARTHQUAKES  (Q) 
AND  PRESUMED  EXPLOSIONS  (EX) 

Case 

No.  of  IF.  *nts 

Parameters 

Probability 
of  Errors 
(percent) 

Q 

EX 

mb 

LTM 

mb-Ms 

1 

it 

12 

X 

B 

X 

B 

6 

2 

it 

12 

B 

X 

5 

3 

ma 

to 

■ 

B 

0 

4 

B 

to 

A 

X 

X 

1 

6 

5 

B 

10 

1 

X 

1 1 

X 

0 

6 

H 

to 

■ 

X 



1 

X 

0 

For  the  first  two  cases,  all  events  of  each  class  were  used.  By  using  only  short-period 
discriminants  the  probability  of  misclassifica+ion  is  ~5  or  6 percent.  The  error  is  reduced 
slightly  by  using  LTM  and  CMPLX  without  m^,  because  m^  is  not  correlated  with  the  other 
parameters  and  simply  adds  noise  to  the  problem.  Cases  3 to  6 utilize  the  data  for  which 
NORSAR  Mg  values  were  available.  As  to  be  expected,  parameters  Mg  and  produce  the 
smallest  probability  of  erior,  essentially  0 in  Case  3.  In  Case  4,  the  short-period  discriminants 
for  the  s^me  event  yield  an  error  probability  of  6 percent.  Finally,  Cases  5 and  6 combine 
(m^-Mg)  with  a single  short-period  discriminant.  Again,  the  overwhelming  effect  of  the  param- 
eter (mj3-Mg)  causes  the  probability  of  error  to  vanisn. 

Because  the  two  classes  of  events  are  snail,  the  statistics  given  in  these  examples  are 
unreliable.  This  is  likely  to  remain  a seric  problem  if  o.ily  short-period  parameters  are 
used.  Strong  variations  in  the  tectonics  and  attenuation  across  Asia  imply  that  data  must  be 
carefully  studied  in  small  regions.  Clearly,  difficulties  in  separation  will  occur  at  lower  m^ 
when  detection  and  noise  problems  become  severe. 

The  best  way  to  stabilize  the  statistics  for  small  classes  is  to  use  several  recording  sites. 
This  mquires  an  adaptation  of  pattern-recognition  techniques  to  multistation  data,  which  is  now 
being  investigated  for  a wide  range  of  magnitudes. 

C.  W.  Frasier 
A.  Shakal 
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E.  DISCRIMINATION  BETWEEN  SOURCE  MECHANISMS 


A study  of  the  feasibility  of  discriminating  among  source  mechanisms  of  mid-ocean  ridge 
eartnquakes  is  under  way.  Since  earthquakes  in  the  mid-ocean  ridge  systems  have  principally 
either  normal  (on  the  ridge  proper)  or  strike-slip  (on  transform  faults)  source  mechanisms,** 

discrimination  between  these  two  types  will  be  attempted  using  the  modified  Karhur.en-Loeve 
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expansion  of  Fukunaga  and  Koontz  as  discussed  above.  Events  of  the  central  Mid -Atlantic 
Ridge  are  being  studied  by  the  analysis  of  single -station  signals,  the  station  being  the  LASA  or 
NORSAR  array. 

Data  from  earthquakes  with  epicenters  near  earthquakes  with  previously  determined  focal 
mechanisms* *~* ^ form  the  observation  vectors  of  the  training  set  for  the  recognition  algorithm. 
The  measurements  making  up  these  observation  vectors  include  those  frequently  used  in 
earthquake-explosion  discrimination  studies,  as  above.  In  addition,  other  discriminants  which 
might  be  more  appropriate  for  this  problem  are  being  considered.  Hopefully,  information 
might  be  obtained  from  this  study  :o  allow  a better  understanding  of  the  wide  variation  seen  in 
typical  Mg-m^  plots  and  to  infer  how  much  of  this  might  be  due  to  real  effects  of  the  mechanism 
of  the  source. 

Certain  tentative  results  suggest  themselves  from  preliminary  analysis.  First,  it  appears 
that  the  separability  of  these  two  populations  is  not  as  good  as  in  the  earthquake-shot  populations. 
Thus,  the  level  of  correct  separation  will  not  be  as  high  at  similar  magnitudes.  Secondly,  it 
appears  that  ridge  earthquakes  have  more  similarity  among  themselves  than  do  transform 
earthquakes.  However,  further  analysis  of  more  data  will  be  necessary  before  solid  conclusions 
can  be  made.  A.  Shakal 
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Fig.  1-1.  Typical  short-period  P-waves 
for  two  classes  of  events.  Class  1 is 
emergent;  Class  2 is  impulsive. 
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Fig.  1-2.  Ms:mb  for  Class  1 (emergent)  events  O and  Class  2 (impulsive) 
events  X.  Multiple  explosion  Ms:  m^,  would  fall  into  box  B from  single 
explosion  region,  box  A.  Its  short-period  record  is,  however.  Class  1. 
Solid  line  is  trend  for  all  Kazakh  earthquakes.^ 
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Fig.  1-3.  Comparison  of  conventional  and  suggested  methods  for  computing  Mg. 
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Fig.  1-4.  Example  of  extraction  of  instantaneous  amplitude  and  frequency  of  surface 
wave  by  means  of  discrete  Hilbert  transform. 
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LOG  THIRD  MOMENT  (LTM) 


QUADRATIC  COMPLEXITT  (CMPLX) 

Fig.  1-5.  Separation  of  presumed  explosions  in  Eastern  Kazakh 
from  nearby  Central  \sian  earthquakes  using  LTM  and  quadratic 
CMPLX. 


Fig.  1-6.  Separation  of  tho  e events 
of  Fig.  1-5  which  have  Measured 
NOP.SAR  Ms  values. 
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II.  FOCAL  DFPTI1 


A.  FOCAL  DEPTH  RESOLUTION  USING  THE  DARPA  SEISMIC  NETWORK 

Preliminary  theoretical  studies  have  been  made  of  the  expected  accuracy  with  which  seismic 
event  locations  will  be  determinable  by  the  DARPA  seismic  network.  Other  calculations  of  this 
same  general  sort  have  been  reported  in  a previous  SATS.*  Estimates  of  the  probable  errors 
contaminating  arrival  time  measurements  are  mapped  into  estimates  of  the  probable  errors  in 
derived  hypocentral  parameters. 

For  the  purpose  of  these  calculations,  stations  which  are  not  yet  installed  were  represented 
by  existing  stations  located  nearby.  Twelve  short-period  stations,  including  both  arrays  and 
seismic  research  observatories  (SROs)  whose  locations  are  known,  were  used: 


LASA 

NORSAR 

SEO 

Seoul,  Korea  (representing  the  Korean  array) 

MSH 

Mashed,  Iran  (representing  the  Iranian  array  and  SRO) 

ANK 

Ankara,  Turkey  (SRO) 

SHL 

Shillong,  India  (SRO) 

WEL 

Wellington,  New  Zealand  (SRO) 

ANP 

Anpu,  Taiwan  (SRO) 

ALQ 

Albuquerque,  New  Mexico  (SRO) 

BOG 

Bogota,  Colombia  (SRO) 

GUA 

Guam,  Mariana  Islands  (SRO) 

AAE 

Addis  Ababba,  Ethiopia  (SRO) 

Figures  II- 1 (a)  and  (b)  show  the  expected  standard  error  of  focal  depth  calculated  on  the 
assumption  that  P-waves  are  detected  by  all  stations  within  95°  of  an  event  and  that  the  arrival 
times  have  standard  errors  of  1.0  sec  (including  the  combined  effects  of  reading  errors,  clock 
errors,  etc.,  and  variations  hetween  actual  travel  times  for  various  paths  and  the  travel  time 
tables  used}.  From  these  figures,  it  appears  that  the  network  may  not  be  able  to  resolve  focal 
depth  well  enough  to  make  it  a useful  discriminant  for  most  shallow  earthquakes  in  Asia. 

B.  R.  Julian 


B.  DETECTING  THE  DEPTH  PHASE  OF  EXPLOSIONS 

Cepstral  analysis,  using  the  Burg  Maximum  Entropy  (BME)  method  to  estimate  the  spectrum 
of  the  log  spectrum  (quefrency).  has  been  used  to  determine  the  arrival  time  of  the  depth  phase 
in  NTS  explosions  recorded  in  Norway  (A  =*  70°).  The  technique  was  described  and  applied  to 
synthetic  data  in  the  last  SATS.2  By  using  the  high-resolution  BME  method,  it  is  possible  to 
maKc  accurate  spectral  estimates  on  the  comp  ex  quefrency  when  the  echo  delay  time  is  small 
and  when  only  that  part  of  the  quefrency  that  contains  significant  power  is  used.  Figures  II-2 
through  II -5  illustrate  the  usefulness  o'  this  technique. 

Figure  II -2  shows  the  short-period  P-wave  of  Boxcar  recorded  at  Oyer  and  the  tapered  ver- 
sion used  in  this  analysis.  The  amplitude  and  phase  components  of  the  spectrum  are  presented 
in  the  top  half  of  Fig.  II -3 . Initially,  the  band  0.3  to  2 sec  was  chosen.  Removing  the  effect  of 
the  instrument  and  attenuation  (t#  = 0.25),  the  linear  trend  in  the  phase  (a  time  delay)  and  the 
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linear  trend  in  the  log  amplitude  (a  source  correction),  and  using  the  window  0.5  to  2 Hz  results 
in  the  complex  quefrency  shown  in  the  bottom  half  of  Fig.  11-3.  Figure  II -4  shows  the  Burg  pre- 
diction error  filter  designed  on  the  complex  quefrency,  and  Fig.  II-5  depicts  the  cepstrum  ob- 
tained from  that  filter.  The  peak  at  0.88  sec  is  interpreted  as  the  delay  time  between  the  P-wave 
and  pP.  Springer*  reports  a pP-P  time  of  0.96  sec  for  Boxcar  using  surface  zero  accelerometer 
records.  The  difference  amounts  to  1 digitizing  umt.  The  delay  time  is  plotted  in  Fig.II-2  on 
the  original  seismogram,  and  the  peak  at  2 sec  in  the  cepstrum  is  also  displayed  as  a time  delay. 
Springer5  reports  a surface  spall  closure  at  3.07  sec  for  'his  event.  The  delay  time  of  2 sec 
found  here  is  consistent  with  the  short-period  waveform  and  may  indicate  that  it  is  not  necessarily 
the  surface  spall  closure  that  produces  teleseismic  energy. 

Similar  results  have  been  found  for  other  NTS  explosions  recorded  in  Norway.  The  shortest 
echo  delay  times  found  so  far  have  been  approximately  0.6  sec.  It  seems  possible  that  this  tech- 
nique will  be  useful  down  to  delay  times  as  short  as  a few  tenths  of  a second  on  data  sampled  at 

20  Hz*  T.  E.  Linders 


C.  CRUSTAL  FOCAL  DEPTH  DETERMINATIONS  FROM  SURFACE  WAVES 

In  the  last  SATS,2  it  was  shown  that  variations  in  upper-mantle  shear  velocity,  particularly 
the  existence  of  thin  high-velocity  layers,  seriously  affect  the  accuracy  of  depth  determinations 
from  Rayleigh-wave  spectra.  This  study  has  been  extended  to  include  the  effects  of  crustal 
structure  variations,  where,  of  course,  thin  high-velocity  layers  may  also  exist.  So  far,  source 
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spectra  have  been  computed  for  three  different  shield  models  - the  Canadian  Shield,  and  Baltic 
and  West  Siberian  Shields,5  where  the  upper-mantle  structure  has  been  assumed  to  be  the  same 
for  all  three.  The  dispersion  curves  for  these  are  identical  to  within  0.05  km/sec  in  the  period 
range  20  sec  upward. 

Spectra  at  a range  of  depths  for  various  fault  orientations  have  been  calculated  wUhin  one 
of  these  models,  and  a search  then  was  conducted  over  depth  for  the  best  fitting  spectra  of  the 
other  models  - a method  identical  with  that  used  in  recent  studies^'7  except  that  the  observa- 
tional spectra  have  been  replaced  by  theoretical  3pectra  for  the  first  model.  Both  phase  and 
amplitude  spectra  have  been  fitted  separately;  for  the  latter,  the  absolute  level  of  the  spectrum, 
governed  by  the  seismic  moment,  has  been  determined  within  the  fitting  scheme  both  from  the 
amplitude  ratio  at  long  periods  and  fi-om  the  average  over  the  whole  period  range  considered 
(16  to  60  sec). 

It  has  been  found,  as  exp  .cted,  that  source  spectra  with  distinctive  features  such  as  "holes" 
or  sharp  phase  changes  can  be  fitted  better  than  those  wUMout  such  discontinuities.  When  the 
source  is  at  depths  too  shallow  for  these  to  appeal  witnin  the  period  range  used,  the  phase  res- 
oluti  - :s  very  poor.  Tests  with  a variety  of  fault  orientations  reveal  no  general  rule  concerning 
the  relative  resolution  or  accuracy  of  phase  and  amplitude  spectra  in  depth  determination.  Fig- 
ure II -6  illustrates  a case  where  the  phase  gives  a much  sharper  fit  than  amplitude;  sometimes 
the  reverse  is  true.  Fitting  phase  and  amplitude  simultaneously  does  not  seem  to  improve  the 
situation. 

Figure  II-7  shows  the  depths  at  which  the  best  fit  was  found  for  spectra  for  the  Baltic  ohield, 
vs  depths  at  which  the  fiited  spectra  were  generated  in  the  Canadian  Shield,  for  a 45°  dip-slip 
fault  fitted  at  azimuths  of  0,  30,  . . . , 150°.  For  other  fault  models,  the  picture  is  generally  the 
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same  - the  effect  of  the  different  earth  models  is  most  serious  below  5 and  above  25  km  depth. 

It  must  be  emphasized  that,  although  this  fit  seems  reasonable,  the  earth  models  used  do  not 
differ  greatly  and  no  errors  in  the  fault  orientation  have  been  considered.  This  study  is  being 
extended  to  include  the  use  of  real  instead  of  synthetic  data. 

R.G.  North 

D.  S RESIDUALS  VS  DEPTH  AT  REGIONAL  DISTANCES  IN  ASIA 

We  have  begun  to  study  the  usefulness  of  the  S-P  interval  in  the  problem  of  depth  determina- 
tion at  near  and  regional  distances.  For  this  study,  we  have  read  the  short-period  recording 
at  World-Wide  Standard  Seismographic  Network  (WWSSN)  stations  of  some  300  Asian  earthquake.;. 
These  events  have  all  been  located  by  the  National  Earthquake  Information  Service  (NEIS)  and 
reported  in  the  Preliminary  Determination  of  Epicenter  (PDE)  lists.  The  chief  difficulty  en- 
countered so  far  is  identifying  the  S arrival,  which  followed  the  same  path  as  the  initial  P ar- 
rival, with  an  accuracy  relevant  to  the  determination  of  depth  within  a few  kilometers  at  shallow 
depth.  A second  problem  lies  in  establishing  reliable  S travel-time  faT,ciions  for  the  region  of 
interest. 

Figure  iI-8  shows  a sample  of  the  results  obtained  to  date.  Here,  we  have  plotted  S travel- 
time anomalies  vs  depth  for  26  events  located  between  7°  and  13°  north  of  New  Delhi,  India  (NDI). 

The  anomalies  (6)  are  related  to  the  S travel-time  function,  T = 10.1  + (A/4.72),  given  for  India 
8 ® 

Ly  Heustis,  et  aL,  and  the  observations  (T  ) by  6 = Tq  - Tc<  The  most  discouraging  aspect  of 
Fig.  II -8  is  that  the  scatter  of  the  data,  over  100  km  for  a given  6,  is  much  toe  great  to  be  of 
use  in  accurate  focal -depth  determination.  Assuming  that  the  NEIS  depths  a*e  accurate,  most 
of  this  scatter  must  be  due  to  our  lack  of  ability  to  consistently  identify  the  initial  S arrival  on 
the  records  read.  A second  disturbing  feature  of  Fig.  II -8  is  that  the  S arrivals  fr<  m near- 
surface  sources  in  a restricted  distance  range  to  the  north  of  NDI  vary  from  4 to  13  sec  later 

g 

than  those  predicted  by  Heustis,  et  aL  The  latter  used  events  located  in  tne  Himalayas  and 
Burma  to  establish  their  S travel -time  curve,  while  in  the  case  of  our  readings  at  NDI  the  events 
were  located  in  the  Hindu  Kush,  Tadzhik,  Kirgiz,  and  western  Sinkiang.  This  indicates  a strong 
path  dependence  on  regional  S velocities  in  Southern  Asia. 

An  encouraging  aspect  of  Fig.  II -8  is  that,  despite  our  lack  of  experience  in  reading  regional, 
short-period  S arrivals,  a definite  trend  with  depth  exists.  This  trend  is  not  a function  of  dis- 
tance since  the  events  are  well  distributed  at  various  depths  in  the  distance  range  studied. 
Whether  we  can  reduce  the  scatter  through  rereading  the  records  and  redefinition  of  the  S travel- 
time function  remains  to  be  shown.  g Dag 

IVLT.  I. in 
J.  R.  Filson 
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Fig.II-4.  Burg  prediction  error  filter 
designed  on  Boxcar's  quefrency. 
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Fig.  II -5 . Cepstrum  for  Boxcar  indicating  that  delay  time  for  pP  is  0.88  sec. 


Fig.  II -6 . Fit  [sum  of  squares  of  difference  between  spectra  at  6 azimuths  (0(30)150°) 
over  period  range  16  to  60  seel  of  spectrum  generated  at  11  km  depth  in  Canadian  Shield 
mor'.el  (CANSD)  for  strike-slip  on  60°  dipping  plane  to  spectra  in  range  0 to  35  km  for 
same  fault  within  Baltic  Shield  (BALTSD).  Minimum  denotes  best  fit.  Depths  indicated 
by  phase  (<f  ),  and  amplitudes  with  moment  determined  by  long-period  spectral  ratio 

[Amp  (Mo^^)|  and  by  averaging  ratio  over  whole  period  range  [Amp  (Mq^^)|  are  13, 

20,  and  15  km,  respectively.  For  this  particular  case,  phase  resolution  is  muchbetter 
than  that  in  amplitude. 


18 


Fig.  11-7 . Best-fitting  depths  within  BAFTSD  for  spectra  generated  in  CANSD  for  45° 
dip-slip  fault  (pure  thrust/normal  fault)  as  function  of  depth.  Depths  given  by  all 
3 fitting  schemes  are  good  except  at  very  shallow  depths  and  below  2 5 km. 


Fig.  11-8.  Short-period  S anomalies 
at  NDI  as  function  of  depth  from  re- 
gional earthquakes  north  of  station. 
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III.  NETWORK  STUDIES 


A.  SEISMIC  DATA  MANAGEMENT  SYSTEM 

Work  is  continuing  on  construction  of  a management  structure  for  seismic  data/information 
resulting  from  the  Integrated  World-Wide  Seismic  System  (IWWSS).  The  extent  and  complexity 
of  the  IWWSS  requires  ti;c  need  of  a documentation  system  capable  of  handling  and  sorting  much 
of  the  detailed  information  and  documents  resulting  from  such  a world-wide  system.  The 
document -handling  utility,  NLS,  currently  in  operation  at  Office-1  on  the  ARPANET,  will  serve 
as  the  means  by  which  text  data  and  information  related  to  the  IWWSS  will  be  sorted,  edited,  and 
stored.  This  facility  at  Office -1  serves  as  the  utility  for  text  editing,  storing,  and  manipulation 
of  the  text  files  that  will  ultimately  become  the  foundation  for  the  data  management  system. 

Once  the  basic  files  a"e  nearly  completed,  ttuy  will  then  be  made  available  to  IWWSS  data  users 
and  related  personnel  via  the  ARPANET  and  the  NLS  system. 

The  acquisition  ol  documents  and  data  related  to  IWWSS  is  dependent  upon  many  of  the  sys- 
tem contractors  and  development  personnel.  Once  information  from  them  can  be  obtained,  it 
will  be  edited  and  incorporated  into  the  data  management  system.  A close  degree  of  cooperation 
among  users  and  developers  will  be  necessary  to  obtain  a system  of  sufficient  detail  to  be  of 
benefit  to  most  users.  Several  organizations  have  been  contacted  and  requested  to  submit  infor- 
mation and  documentation  for  inclusion  into  the  system.  Documentation  that  currently  exists 
for  IWWSS  components  which  have  reached  some  decree  of  stability  are  to  be  included  in  the 
management  system.  One  major  problem  at  present  is  the  la  '•ge  degree  of  fluidity  of  many  of 
the  current  plans  for  IWWSS.  Because  of  this  condition,  it  is  not  possible  to  include  many  of 
the  major  components  of  IWWSS  at  this  time.  Some  components  that  are  firmly  established, 
such  as  datacomputer  datalanguage,  hnve  already  been  included  into  the  management  system. 

The  document  "Datacomputer  Version  0/10  User  Manual"  has  been  entered  in  the  NLS  system 
and  is  available  for  use  within  NLS.  Sorr  restructuring  of  this  manual  was  necessary  to  make 
full  use  of  the  capabilities  of  NLS.  The  document  index,  given  in  the  first  branch  cf  the  file, 
provides  rapid  access  to  any  major  section  i ' that  user  manual.  Once  the  user  determines  the 
appropriate  chapter,  section,  or  topic  of  interest  from  the  index,  the  Command  "Print  Branch" 
or  "Jump  (to)  Link"  will  access  that  area  of  the  document.  The  index  of  this  file  together  with 
those  of  all  other  files  are  contained  in  a master  index  file,  and  will  be  presented  to  the  user 
when  he  enters  the  NLS  system. 

Users  of  the  data  management  system  will  not  necessarily  be  well  versed  in  NLS  command 
language.  To  make  this  information  system  usable  by  the  endre  community,  it  will  be  necessary 
to  ask  the  user  (when  he  enters  NLS)  if  he  is  familiar  with  NJJ5.  If  he  is  not  acquainted  with  the 
NLS  language,  then  he  will  be  presented  with  a brief  introduction  to  the  NLS  and  IWWSS  systems 
and  a subset  of  NLS  commands  needed  for  a "rend  only"  knowledge  of  the  system.  Also,  in  order 
to  make  the  NLS  system  seem  "cleaner,"  i.e.,  less  cluttered  by  heralds,  prompts,  and  echoes, 
these  features  are  eliminated  for  the  new  user.  It  now  appears  that  a user  knowledge  of  only 
8 to  10  commands  will  be  sufficient  to  acquire  a "read  only"  cap  ability  in  NLS. 

Several  files  have  been  set  up  and  now  contain  the  basic  information  and  categories  of  IWWSS 
topics.  These  will  be  the  foundation  of  the  management  system,  and  will  be  expanded  and  devel- 
oped as  more  information  becomes  avail aole.  In  addition  to  the  index  file  there  is  an  informa- 
tion file  that  contains  a general  description  of  each  file  in  the  management  system.  Other 
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currently  available  files  include  one  with  the  personnel  names,  organizations,  and  organizational 
responsibilities  for  the  IWWSS.  One  other  file  is  dedicated  to  a description  of  all  SROs  and  all 
other  seismic  sites  that  make  up  part  of  the  IWWSS.  A complete  site  description  and  location 
will  be  provided  along  with  descriptions  of  the  instrumentation  at  each  site  as  soon  as  these  sites 


become  operational. 


R.M.  Sheppard 


B.  DETECTION  CAPABILITY  OF  A 32 -STATION  SUBSET 

OF  ISM  STATIONS  DURING  THE  ISM  PERIOD 

An  investigation  into  the  detection  capability  of  a 32-station  subset  (Fig.III-1)  of  the  approx- 
imately 150  stations  used  by  the  ISM  experiment*  has  been  completed.  A simple  statistical 
model  to  predict  the  probability  of  an  N-station  detection  from  an  M -station  network  was  devel- 
oped. The  data  from  the  32-station  set  were  used  to  estimate  the  single-station  detection  prob- 
ability for  events  with  4.6  < m^  <4.7  and  withm  90®  of  the  station.  This  probability  was  entered 
into  the  model,  and  a network  detection  probability  table  was  generated.  A 1 echnical  Note  re- 
porting the  detailed  results  of  this  investigation  is  in  preparation. 

The  selection  of  the  32-station  subset  was  ad  hoc  to  give  good  global  coverage  while  using 
stations  which  contributed  data  to  a large  number  of  ISM  events.  The  ISM  event  list  is  composed 
of  996  events  W'hich  occurred  between  20  February  and  19  March  1972.  The  subset  stations  de- 
tected and  located  819  of  these  events.  Many  of  the  177  *■"  ents,  reported  on  the  ISM  list  but  not 
reported  on  the  subset  list,  were  small  events  requiring  a more  compact  geographically  distrib- 
uted network  than  the  subset  network,  to  meet  the  epicenter  location  criteria  used  by  this 
experiment. 

A comparison  was  made  between  the  USC.S  PDF  list  and  the  event  list  created  using  the 
32 -station  subset  network.  The  subset  list  showed  508  events  not  reported  by  the  PDE  list,  while 
the  PDF'  list  reported  43  events  not  included  in  the  subset  list.  All  the  stations  used  in  the  com- 
putation of  events  listed  by  the  subset  network,  with  the  exception  of  CIIG  and  YKA,  were  contrib- 
utors to  the  PDE  computations.  However,  the  detection  and  reporting  level  to  USGS  and  for  our 
experiment  were  quite  different.  For  example,  the  detection  and  reporting  level  for  the  USGS 
and  for  the  ISM  subset  were  compared  for  stations  MBC  and  KRL.  The  additional  detections 
available  to  Lincoln  vastly  improved  the  stations'  detection  capability  for  events  with  m^  < 4.5. 

Using  the  32-station  subset  network  of  well  read  and  reported  stations,  a single-station  de- 
tection probability  of  0.63  was  estimated  for  events  with  4.6  < mb  < 4.7  and  90°.  Table  III- 1 
shows  how  this  converts  into  network  detection  probability.  For  example,  if  20  stations  of  the 
subset  are  within  90°  of  an  event  m^  of  4.6,  the  probability  is  100  percent  that  5 of  these  stations 
will  detect  it  and  the  probability  is  greater  than  90  percent  that  10  stations  of  the  subset  will  de- 

tec*  ***  R.E.  Needham 

C.  SEISMIC  EVENT  LOCATION  PROCEDURES 

The  event-locatiwn  program  developed  for  the  ISM  project^  is  undergoing  continuing  modifica- 
tion and  improvement.  These  changes  include:  (1)  the  addition  of  new'  seismic  phases,  (2)  the 
ability  to  restrain  latitude  a-id/or  longitude  and/or  origin  time,  as  well  as  focal  depth,  and 
(3)  provision  for  atxempting  to  automatically  identify  and  reject  observations  inconsistent  with 
the  other  data.  This  last  provision  is  particularly  useful  because  data  are  frequently  contam- 
inated by  a few  gross  errors,  the  detection  of  which  consumes  i great  deal  of  analysts'  time. 
Tables  III-2(a)and  (b)  show  an  example  from  the  ISM.  an  earthquake  in  the  southern  Kuril  Islands. 
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TABLE  ni-2(a) 
NITIAL  (ISM)  LOCATION 
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TABLE  m-2(b) 
FINAL  LOCATION 
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Although  there  is  nothing  obviously  wrong  with  the  epicenter  determined  for  the  ISM,  a much 
better  fit  to  the  data  (and  an  18-km  shift  of  the  epicenter)  is  obtained  if  the  arrival  time  data  from 
NORSAR  are  discarded. 

A Technical  Note  describing  the  location  program  in  detail  is  in  preparation. 

B.  R.  Julian 


D.  STABILITY  OF  TELESEIi'MIC  TRAVEL-TIME  RESIDUALS 

It  has  been  shown  that  seisi.  ic  wave  velocities  in  seismic  regions  may  change  with  time, 

and  one  way  of  measuring  tl  is  vacation,  particularly  when  a dense  local  netwrork  is  unavailable, 

is  to  study  teleseismic  P-wave  tnve.-time  residuals  at  a particular  station  to  see  whether  these 

1 4 

are  subject  to  temporal  variation.  ' Such  changes  in  travel  time  will  obviously  have  a measur- 
able effect  on  the  location  ability  of  a seismic  network.  Considerable  care  must  be  exercised 
to  ensure  that  the  variations  observed  (~0.2  to  0.4  sec)  do  not  arise  from  epicenter  migration, 
location  errors,  station  timing  errors,  and  many  other  factors,  and  in  general  only  deep  events 
are  used  and  the  residual  averaged  over  intervals  of  about  6 months.  A suitable  way  to  test  the 
stability,  or  lack  of  it,  of  residuals  is  to  study  those  determined  at  teleseismic  distances  for 
nuclear  explosions  at  the  Nevada  Test  Site  (NTS).  Here  the  location,  origin  time,  and  depth  are 
known  for  over  100  useful  events  over  the  last  10  years,  and  the  only  errors  left  are  those  in 
station  timing  and  reading  and  possibly  from  movement  of  shot  location  within  NTS. 

For  20  stations  which  have  regularly  reported  arrivals  from  NTS  events  over  1964  to  1972, 
the  travel-time  residuals  calculated  by  the  International  Seismological  Center  (ISC)  with  respect 
to  the  locations  given  by  the  AEC  were  plotted  as  a function  of  time.  An  immediately  noticeable 
feature  of  the  results  W'as  a variation  from  station-to- station  of  the  spread  and  standard  devia- 
tion of  the  residuals.  This  could  be  due  to  reading  or  timing  errors,  grossly  heterogeneous  sub- 
station structure,  or  to  "kinks"  in  the  travel-time  curve.  Stations  giving  consistent  residuals 
with  a standard  error  of  less  than  0.3  sec  included  UME,  SOD,  ARE,  GIL,  NP-  and  KTG. 
Standard  errors  in  excess  of  0.5  sec  were  observed  for  TRO,  KON,  GRR,  FLN,  SSF,  and  COL 
(Fig.III-2). 

Larger  residuals  (earlier  arrivals)  were  almost  invariably  observed  for  the  larger  shots, 
either  because  of  better  time  picks  or  the  grea.er  shot  depths  used.  Shot  elevations  (above  sea 
level)  range  from  0.5  to  1.9  km,  introducing  a non-negligible  travel-time  correction  which  was 
not  included  by  the  ISC.  By  assuming  a near-surface  velocity  of  2.5  km/sec  to  correct  for  depth, 
these  times  could  be  moved  closer  to  the  mean.  A study  of  the  correlation  of  location  within 
NTS  with  residual  values  yielded  negative  results. 

At  two  stations,  ADK  (Adak)and  MAT  (Matsushiro),  apparent  temporal  variations  in  residual 
values  were  observed  both  of  ~0. 5-sec  magnitude  and  durations  2 to  3 years.  These  are  shown 
in  Fig.III-3,  the  variation  in  both  cases  commencing  in  late  1967.  A comparison  of  residuals 
at  pairs  of  close  stations,  such  as  (GRR-FLN),  (FLN-LOR)  (France)  reveals  that  the  standard 
error  of  the  difference  in  residual  between  them  is  significantly  less  than  the  error  in  that  for 
the  individual  stations.  This  indicates  that  heterogeneity  near  some  stations  may  be  the  cause 
of  the  large  standard  errors  observed.  In  some  instances,  e.g.,  COL,  this  i - obviously  not  the 
case,  since  at  a ver^'  close  station  (GIL)  the  standard  deviation  of  the  residual  was  only  0.24  sec 
compared  with  0.54  sec  for  the  former.  This  can  only  be  ascribed  to  bad  time  picks,  and  em- 
phasizes that  considerable  caution  must  be  exercised  in  using  bulletin  data. 

R.G.  North 
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Fig.  III-l.  Locations  of  stations  used  in  32-statiot.  network. 
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Fig.  III-2.  F-wave  travel-time  residuals  from  NTS  at  ARE  (Arequipa,  Pe  m, 
A = 67.8")  and  COL  (College,  Alaska,  A = 33.5°).  ARE  readings  appear,  wuh 
standard  deviation  of  0.25  sec,  to  be  much  better  than  those  at  COL  (a  - 0.54). 


Fig.III-3.  Apparent  time  variation  of  residuals  fro  n NTS  at  ADK  (Adak,  Aleutians, 
A -44.4°)  and  MAT  (Matsushiro,  Japan,  A = 79.6#).  In  both  cases,  residual  appears 
to  deviate  b\  about  0.5  sec  from  mean  for  2-year  period  commencing  in  mjd-l°67. 
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IV.  EARTH  HETEROGENEITY 


A.  EVIDENCE  FOR  ANISOTROPIC  SCATTERING  IN  THE  UPPER  MANTLE, 

DEDUCED  FROM  P,  pP,  PcP,  AND  pPcP  PHASES 

In  a previous  study  of  the  core  mantle  boundary,  we  observed  that  short-period  PcP/P 
amplitude  ratios  measured  at  LASA  are  a strong  function  of  focus  depth.  Shallow  events  from 
0-  to  100-km  depth  have  PcP/P  ratios  3 to  4 times  the  ratios  for  events  deeper  than  200  kir . 

Our  conclusion  was  that  a strong  scattering  mechanism  in  the  upper  mantle  rather  than  inelastic 
attenuation  was  producing  the  observed  depth  effect  of  PeP/P  ratio.  This  scattering  is  not  i -jo- 
tropic  but  is  a strong  fur^Uon  of  the  horizontal  phase  velocity  of  the  rays  passing  through  the 
upper  mantle.  Rays  of  rU^h  dT/dA  are  attenuated  more  than  waves  of  low  dT/dA  along  compa- 
rable path  lengths. 

Additional  data  nave  been  obtained  which  clearly  isolate  the  scattering  effect  to  the  upper 
mantle  above  the  focus.  This  s done  by  using  the  surface-reflected  phases  pP  and  pPcP  in 
conjunction  with  P and  PcP  phases.  Figure  IV-1  shows  the  ray  paths  for  these  phases  for  a 
particular  distance  and  focus  depth.  In  order  to  isolate  the  upper-mantle  ray  paths,  we  computed 
the  ratio  of  pPcP/pP  to  PcP/P  amplitudes  for  20  LASA  events.  There  are  several  advantages 
to  using  this  ratio: 

(1)  The  takeoff  angles  of  pPcP  and  pP  are  very  close  at  teleseismic  distances,  and  one 
may  assume  that  similar  source  amplitudes  are  radiated  along  each  direction.  Sim- 
ilarly, one  can  iume  equal  source  radiation  along  the  ray  paths  of  PcP-  and 
P-waves.  By  taking  (pPcP/pP)/(PeP/P)  ratios,  we  thus  eliminate  the  unequal 
source  radiation  in  the  up-and-down  directions  from  the  focus. 

(2)  Reflection  coefficients  at  the  free  surface  divide  out  because  the  incident  angles 
of  pPcP  and  pP  are  nearly  equal. 

(3)  pPcP  and  PcP  rays  sample  the  same  part  of  the  core  mantle  boundary  and  experi- 
ence the  same  spherical  spreading  as  a function  of  distance  and  depth.  Also,  pP 
and  P sample  the  same  part  of  the  lower  mantle.  Thus,  the  ratio  (pPcP/pP)/ 
(PcP/P)  cancels  the  spherical  spreading,  attenuation,  and  core  mantle  boundary 
reflection.  This  essentially  eliminates  the  earth  structure  below  the  source  and 
isolates  the  effects  of  the  mantle  above  the  source. 

Figure  IV-2  shows  the  (pPc P/pP)/( PcP/P)  amplitude  ratios  for  20  events  recorded  at  LASA 
as  a function  of  epicentral  distance.  The  amplitude  ratios  vary  by  almost  an  order  of  magnitude 
with  no  trend  with  distance.  The  same  data  are  plotted  as  a function  of  focus  depth  in  Fig.  IV-3. 
Except  for  one  data  point,  there  is  a clear  trend  of  this  ratio  from  1 for  the  shallowest  event, 
to  3 or  greater  at  depths  >200  km.  The  values  are  scattered  but  significantly  >1  in  all  cases 
except  one.  As  previously  discussed,  these  amplitude  ratios  should  cancel  out  the  earth  model 
effects  below  the  source,  and  isolate  the  ray  paths  above  the  source. 

Figure  IV -4  shows  the  theoretical  effect  of  Kanamori's  Q model2  on  (pPeP/pP)/(PeP/P) 
ratios  as  a function  of  depth,  distance,  and  peak  frequency  of  the  phases.  At  1 Hz.,  the  effect  of 
attenuation  is  to  elevate  the  ratio  smoothly  with  depth  to  values  >1.  This  effect  increases  with 
depth  of  focus,  depending  on  the  epicentral  distance  as  shown  by  the  solid  lines  for  40°  and  80°. 

At  0 Hz,,  no  attenuation  would  occur  and  the  region  would  collapse  into  a single  line  of  value  1 
fo'’  aU  depths  and  distances. 


The  effect  of  attenuation  was  removed  from  the  data  using  Kanamori's  Q model.  The  cor- 
rected (pPcP/pP)/(pP/P)  ratios  are  shown  in  Fig.  IV-5.  If  the  uncorrected  ratios  were  explained 
by  such  a low  Q model,  then  the  correeted  ratios  in  this  figure  would  scatter  about  1.0  for  all 
depths.  However,  the  corrected  ratios  are  still  significantly  above  1.  In  order  to  correct  the 
ratios  to  a level  of  - 1 for  all  depths,  an  average  Q of  <40  for  the  upper  400  km  of  the  mantle 
would  be  required. 

Such  low  Q values  for  short-period  data  are  clearly  unacceptable,  une  explanation  is  an 
anisotropic  scattering  mechanism,  which  attenuates  the  ray  paths  above  the  source  for  pP  phases 
more  severely  than  pPcP  phases.  As  previously  suggested,1  thin  horizontal  laminations  of  low 
rigidity  sandwiched  between  thick  layers  of  solid  material  could  produce  such  a phenomenon. 


Theoretical  studies  of  scattering  of  P-waves  through  more  realistic  mantle  models  are  being 


undertaken. 


C.  W.  Frasier 


D.  K.  Chowdhury 


B.  SCATTERING  OF  ELASTIC  WAVES  IN  A MEDIUM  WITH  ANISOTROPICALLY 
DISTRIBUTED  RANDOM  INHOMOGENEITIES 

Previous  theoretical  studies  of  elastic  wave  scattering^  ^ in  random  media  have  dealt  only 
with  the  case  in  which  the  inhomogeneities  are  distributed  without  regard  to  spatial  orientation. 
However,  recent  studies  of  PcP/P  amplitude  ratios1, ^ suggest  that  in  the  upper  mantle,  P-wave 
scattering  depends  strongly  on  the  angle  of  incidence,  an  effect  which  cannot  be  explained  by 
these  isotropic  theories.  The  results  of  these  studies  also  indicate  that  a significant  fraction 
of  the  energy  (50  percent  or  more)  in  short-period  P—vaves  may  be  scattered  in  a single  passage 
through  the  upper  mantle,  a phenomenon  potentially  capable  of  seriously  biasing  many  types  of 
seismic  investigations. 

4 5 

We  have  extended  the  theory  given  by  Knopoff  and  Hudson  ’ to  the  case  of  a medium  with 
an  anisotropic  distribution  of  random  inhomogeneities.  That  is,  the  elastic  constants  X and  p. 
and  the  density  p are  assumed  to  be  functions  of  position  r given  by 

X = \Q  + 6X(  r ) ; p = pQ  + <5p(  r ) ; p = pQ  + 6p  ( r ) 

where  6X,  6p,  and  <5p  are  specified  only  in  terms  of  their  correlation  function: 

NX(  rr  r2)  = <6X(  r1),  fiXfrM)  NXp(  7^ . rvj  = <flX(  ),  <5p(r2)> 

Np(  Fj , r^)  = <6p(~1 ),  op(r2)>  Npp ( , F, ) = <fip(  r^  ),  6p(r^)} 

NpCFi,  r2)  = <6p(  ?i ).  <5p(r2)>  NXpfr, . r2)  = <6X(  ).  6p  (Tv^  )> 

where  the  brackets  indicate  averaging  over  the  ensemble  of  all  possible  inhomogeneities. 

4 5 

Previous  studies  such  as  Kr.opoff  and  Hudson  ' have  assumed  these  correlations  to  be  functions 
of  |r^  - r2|,  usually  of  Gaussian  form: 

N(7^,  r^)  = a2  exp 
whereas  we  have  used 
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Figure  IV -6  gives  some  preliminary  results  for  P-to-F  scattering,  comparing  the  isotropic 
case  with  the  case  in  which  a = b = 5c.  It  is  evident  from  these  results  that  anisotropy  in  the 
distribution  of  inhomoger.oities  profoundly  effects  the  scattering  of  elastic  waves,  and  in  fact 
that  the  most  energy  is  scattered  in  a direction  controlled  by  the  "grain"  of  the  medium,  rather 
than  by  the  direction  of  the  incident  wave. 

A paper  containing  more  details  of  the  theory  and  compiling  results  for  a variety  of  cases 
(including  P-to-S,  S-to-P,  and  S-to-S  scattering)  is  in  preparation. 

W.  B.  Lee 
B.  R.  Julian 


C.  SCATTERING  IN  THE  MID-TOUPPER  MANTLE 

Both  forward  and  backward  scattering  by  random  inhomogeneities  within  the  crust  have  been 

suggested  as  sources  of  short-period  seismic  codas/*'7  Forward  scattering  at  or  near  the  core- 

8 9 

mantle  boundary  has  been  invoked  as  a soui  ce  of  precursors  to  PKKP  and  PKP.  ’ Here  a simple 
model  for  teleseismic  scattering  is  suggested,  and  evidence  is  presented  that  strong  scattering 
occurs  in  the  mid-to-upper  mantle. 

Since  multiple  scattering  greatly  reduces  the  scattered  energy  in  a given  direction,  a model 
is  assumed  in  which  singly  scattered  waves  contain  the  major  scattering  energy.  To  isolate 
scattered  energy  from  mantle  regions,  short-period  codas  between  major  phases  from  shadow 
zone  events  have  been  analyzed.  This  technique  greatly  reduces  the  complexity  of  the  seismo- 
grams to  be  studied.  Figure  IV -7  represents  a cross  section  of  the  earth  showing  the  shadow 
zones  for  two  events  that  are  more  than  95°  apart,  at  points  A and  B.  The  ray  AP^  is  the 
P-wave  path  that  grazes  the  core  and  defines  the  strictly  geometrical  core  shadow.  The  ray 
AP^  defines  the  last  observable  ray  path  of  P diffracted  in  the  geometrical  shadow  zone  for  a 
source  at  A and  a receiver  at  P^.  BP^  and  represent  the  same  paths  for  a source  at 

point  B. 

Assuming  a receiver  at  point  A is  used  to  study  an  event  fr^m  point  B,  it  is  evident  that 
only  inhomogeneities  that  lie  within  the  intersection  oi  the  regions  above  the  two  true  shallow 
zones  can  contribute  singly  scattered  signals.  The  arrival  time  of  scattered  waves  at  point  A 
for  each  scattering  point  (e.g.,  points  1 and  2 in  Fig.  1V-7)  may  be  calculated  as  the  sum  of  the 
wavefront  times  which  intersect  at  that  point.  This  arrival  time  plus  the  observed  slowness  at 
A define  the  location  of  the  seatterer.  For  example,  scattering  at  points  1 and  2 would  result 
in  scattered  travel  times  of  17  minutes  for  the  configuration  of  source  and  receiver  shown,  hut 
waves  from  these  two  points  would  arrive  at  point  A with  quite  differ  t slownesses. 

Of  course,  the  model  is  not  restricted  to  only  P-to-P  scattering;  P-to-S  and  S-to-P  scat- 
tering might  also  occur.  Referring  to  [joints  1 and  2 again,  and  assuming  a v^/v^  ratio  of  2, 
the  travel  time  to  point  A of  the  scattered  waves  would  be  27  minutes. 

Figure  IV-8  shows  a possible  application  of  this  model  to  results  obtained  from  the  .hady 
of  the  two  following  shallow  events  near  the  Philippine  Islan  is: 

Azimuth 

i 05  ° 305° 

1 07  0 303° 
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Short-period  data  from  the  lASA  E-  and  F-rings  were  bandpass  filtered  with  a window  of 
0.9  ± 0.4  Hz.  These  data  were  then  used  to  form  Vespagrams.  In  both  cases,  an  arrival  was 
found  at  17  minute.'.  ± 4 sec  after  the  event  time  with  a slowness  of  6.4  sec -deg  . This  slowness 
corresponds  to  a P-wave  arrival  from  a distance  of  64°.  The  intersection  ol  this  ray  with  the 
17 -minute  isotime  line  defines  the  position  of  the  proposed  scatterer.  This  places  the 
scatterer  at  a depth  of  about  1000  km  under  the  intersection  of  the  American  and  Pacific  plates 
and  the  great  circle  path.  Analysis  of  a third  event  from  the  same  region  but  at  a greatei  dis- 
tance (121°)  also  provided  an  arrival  appropriate  to  this  scatterer.  Or  course,  U.e  energy  in  this 
arrival  could  be  equally  well  due  to  multiple  scattering  in  the  receiver  region  of  a late  P-phase. 
This  will  be  resolved  by  undertaking  frequency -wavenumber  spectral  analysis  of  these  events. 

J.  Scheimer 

A,  V.  Nikolaev 

T.  E.  Landers 

D.  SHEAR  WAVE  VELOCITY  IN  THE  EARTH'S  LOWER  MANTLE 

Unlike  compressicnal  velocity,  the  shear  velocity  in  the  earth's  mantle  is  known  so  far  at 
best  within  an  uncertainty  of  0.1  km/sec  or  more.  Recent  observations  of  large  numbers  of 
overtone,  of  normal  modes  have  given  us,  however,  a better  estimate  of  shear  velocity  (sec, 
for  example,  Gilbert  and  Dziewonski*  °).  Precise  estimates  for  radial  variation  in  shear  velocity 
have  become  important  now  more  than  ever  as  present  knowledge  of  the  earth  is  getting  more 
complex  radially  and  laterally.  Any  kind  oi  thermally  activated  inhomogeneity  would  affect  shear 
velocity  more  than  comprcssional  velocity,  and  thus  shea"  velocity  variation  would  be  a sensitive 
indica..or  for  the  presence  of  thermal  p’umes  or  mantle-wide  thermal  convection.  At  the  present 
time,  body  wave  travel  time  and  amplitude  data  are  more  effective  tools  for  studying  lateral 
.leterogenoity  in  the  earth’s  mantle  than  the  normal  mode  data.  Therefore,  we  need  a body  wave 
standard  tor  shear  velocity,  eliminating  any  possibility  of  inconsistency  among  normal  mode  data 
and  body  w ve  data  due  to  diiferer.e  e in  the  wavelength  irvolved. 

A standard  for  radial  variation  in  shear  velocity  has  been  obtained  from  the  inversion  of 
deep  shock  travel-time  data,  thus  eliminating  the  effect  of  near-source  heterogeneity  in  the  upper 
mantle  and  enhancing  the  reading  accuracy.  Our  velocity  model  is  presented  in  Fig.  IV -9  and  is 
tabrlated  in  Table  IV-1.  Nonuniqueness  in  this  velocity  model  has  been  studied  in  the  framework 
Oi  the  inversion  theory  of  Backus  and  Gilbert.**  According  to  this  theory,  the  vertical  error  bar 
(Fig.  I V - 9 ) in  our  velocity  model  (termed  "spread")  results  from  the  finiteness  in  the  number  of 
the  observed  data  and  gives  a measuia  of  esolving  length.  The  horizontal  bar  (Fig.  1V-9),  on 
the  other  hand,  arises  due  to  error  in  the  observation  and  also  due  to  the  inadequacy  of  the  model 
in  fitting  the  observed  data.  These  two  errors  - spread  and  error  in  velocity  - : re  not  independent 
but  one  can  be  reduced  at  the  expense  of  the  other.  The  functional  relationship  between  the  two 
errors,  commonly  known  as  "trade-off  curves,"  is  shown  in  F’g.  lV-10  for  different  depths  in 
the  lower  mantle.  One  immediate  application  of  these  trade-off  cur  ves  is  to  assert  that,  with  our 
data,  one  cannot  parameterize  the  shear  velocity  model  with  50-km -thick  and  resolvable  layers 
ii.  the  mantle.  Such  a parameterization  would  lead  to  an  error  in  the  computed  travel  times  of 
a larger  magnitude  than  observed.  The  resolving  kernels  for  corresponding  depths  are  shown 
in  Fig.  IV -11.  These  resolving  kernels  show  that  error  in  the  assumed  shear  velocity  distribution 
in  the  upper  mantle  would  not  affect  the  velocity  model  in  the  lower  mantle  below  1000  km.  Apart 
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from  the  reliable  estimates  of  velocity,  other  features  of  our  shear  velocity  model  include  the 
absence  of  velocity  inversion  near  the  core  mantle  boundary,  and  fluctuations  of  long  wave- 
length throughout  the  mantle.  Among  all  the  ex  sling  shear  velocity  models,  our  model 
comes  closest,  to  the  i066  B model  of  Gilbert  and  Dzicwonsld*  ^ which  was  obtain°d  from  1064 


normal  mode  fundamental  and  overtone  eigenperiods. 


M.  K.  Sengupta 
B.  R.  Julian 
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Fig.  IV -1.  Ray  paths  for  P,  pP,  FcP, 
and  pPcP  phases  recorded  at  same  epi- 
centra)  distance.  L.V.Z.  denotes  low 
velocity  zone  of  upper  mantle. 
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DISTANCE  PROM  LAS  A (deq) 

Fig.IV-2.  Amplitude  ratios  (pPc P/'pP)/(PcP/P)  for  20  events  recorded 
at  LASA.  Open  circles  denote  earthquakes  south  of  LASA,  wherear  black 
circles  are  events  from  north. 
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Fig.  IV-4.  Effect  of  attenuation  on  (pPcP/pP)/(PcP/P)  ratios  as  a function 
of  distance,  deptli  of  focus,  and  frequency.  Kanamori's  Q model  (1967)  is  used. 
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Fig.  IV-5.  Amplitude  ratios  (pPeP/pP)/( Pep,  P’*  of  Fig.IV-2  corrected 
for  Kanamori's  Q model. 
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Fig.  IV -6.  Directional  dependence  of  amplitude  of  P-waves  scattered  f rom  inc  ident 
P-wavc  in  case  where  correlation  lengths  are  5 times  greater  in  horizontal  (x  and  y) 
directions  than  in  vertical  (z)  direction.  Arrows  indicate  directions  of  incident 
waves  (0°,  30°,  60°,  and  90°  to  +x  axis).  B'or  comparison,  dashed  line  gives  cor- 
responding result  for  isotropic  case. 


Fig.  IV-7.  Shadow  zones  for  two  events  >95°  apart.  Rays  P.  and  P define 

/V  fc>  ( j 

strictly  geometrical  shadow  zones  for  event  s at  points  A and  13.  Rays  P ^ 
and  P^j  are  typical  paths  for  diffracted  P-waves  from  these  sources.  Arcs 
represent  position  of  P-wave  wavefront  from  each  source  at  successive 
1 -minute  intervals.  Thus,  t^ , t^,  t^  , and  t^  ^ define  positions  of  P-wave 
at  1,  2,  *,  and  10  minutes,  respectively,  after  event  at  point  A. 


Fig.  IV-8.  Model  for  earthquake  discussed  in  text.  Line  LL  is  "iso-time"  li;u 
which  represents  locus  of  possible  positions  of  P-P  seatterers  that  will  allow 
total  travel  time  between  points  A and  B.  Intersection  of  this  locus  with  ray 
defined  at  LASA  array  by  slowness  of  6.4  gives  probable  position  of  hypothesized 
single  scatterer  in  mantle. 


Fig.  IV-9.  Snear  velocity  model  for  earth's  lower  mantle.  Horizontal 
and  vertical  error  bars  are  explained  in  text. 


V.  SOURCE  MECHANISM  STUDIES 


A.  P-WAVEFORMS  FROM  SHALLOW  EARTHQUAKES 

A significant  advance  might  bo  made  in  discrimination  seismology  if  we  could  accurately 
interpret  the  P-waveform,  both  long  and  short  period,  from  shallow  earthquake  sources.  Here 
"long  and  short  period"  refers  to  the  response  of  the  typical  World-Wide  Standard  Seismographic 
Network  (WWSSN)  instruments.  The  interpretation  of  the  short-period  seismograms  is  more 
relevant  to  the  discrimination  problem  because  of  their  higher  magnification;  however,  the 
long-period  data,  because  of  the  longer  wavelengths  involved,  are  less  contaminated  by  minor 
variations  in  structure  along  the  path.  We  propose  to  use  the  long-period  data  as  a crutch  to 
help  us  identify  the  gross  features  of  the  source  and  to  help  us  interpret  salient  features  on  the 
short-period  records.  In  Fig.  V-l(a),  we  demonstrate  that  this  is  a valid  way  to  proceed. 

Trace  1 is  the  long-period  vertical  recording  of  a snallow  Asian  earthquake  recorded  at 
Berkeley,  California  (BKS).  This  trace  was  deconvolved  witn  the  proper  instrument  response 
to  estimate  ground  displacement,  and  then  reconvolved  with  the  short-period  impulse  response 
to  yield  trace  2.  Trace  3 is  the  actual  short-period  vertical  recording  at  BKS  of  this  event; 
most  of  the  salient  features  of  this  trace  are  reproducible  from  the  long-period  record.  The 
assertion  that  relevant  source  information  exists  in  short-period  data  is  supported  in  Fig.  V-t(b) 
where  the  short-period  traces  from  this  event,  recorded  at  two  European  stations  (NUR  and 
STU),  are  compared.  These  stations  are  at  about  the  same  azimuth  from  the  source,  but  vary 
in  epicentral  distance  by  15°.  The  first  few  cycles  of  these  two  waveforms  are  quite  similar, 
which  indicates  that  their  basic  form  is  due  to  the  source,  but  they  are  dissimilar  from  the  BKS 
short-period  data  recorded  at  a different  azimuth.  The  long-period  P-wave  recordings  of  this 
event  showed  strong  but  regular  vai  .ations  as  a function  of  azimuth. 

Synthetic  P-wave  seismograms  have  been  computed  by  many  workers  (notably  Douglas, 
et  al.,1  and  Helinberger2),  and  examples  of  these  computations  are  shown  in  Fig.  V-2(a-c). 

The  upper  trace  in  Fig.  V-2(a)  is  the  far-field  P-wave  radiation  due  to  a circular  strike-slip 
fault  computed  from  a moJel  developed  by  Savaged  The  traces  beneath  ihis  pulse  represent  it 
on  short-  and  long-period  WWSSN  seismometers.  In  Fig.  V-2(b),  the  middle  trace  represents 
radiation  from  the  same  fault  with  a surface  reflection  added,  and  the  upper  trace  in  (c)  repre- 
sents a multiple  event  of  the  same  type  without  surface  reflections.  Beneath  each  of  these 
source  functions  are  the  attendant  short-  and  long-period  seismograms.  The  point  here  is  that 
one  can  qualitatively  match  the  general  features  of  the  data  shown  in  Fig.  V-l(a-b)  by  the  com- 
bination of  proper  surface  reflections  and,  in  some  cases,  multiple  sources. 

Although  I lon't  feel  confident  enough  in  the  model  used  to  give  a more  specific  interpre- 
tation at  this  time,  I can  make  three  general  observations:  (1)  The  data  support  earthquake 
models  in  which  rupture  begins  at  a single  point  rather  than  simultaneously  along  a line,  (2) 
ambiguity  will  exist  in  most  interpretations  of  source  detail  based  on  one  or  a few  stations, 
and  (3)  earthquakes  with  discrete  multiple  sources  are  common. 

J.  R.  Filson 

B.  POLARIZATION  OF  LONG- PERIOD  BODY  WAVES  FROM  CANNIKIN 

4 5 

A three -component  polarization  program  written  by  Madariaga  ’ has  been  modified  for  the 
data-analysis  console  at  Linco'n  Laboratory,  and  is  being  used  to  investigate  the  polarization  of 
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waves  within  codas  of  earthquakes  and  explosions.  The  method  avoids  the  frequency  domain 
calculations  of  the  polarization  scheme  of  Shimshoni  and  Smith/'  and  is  isier  to  interpret  than 
the  REMODE7  program.  Basically,  this  method  calculates  particle  motion  ellipses  as  a func- 
tion of  time  by  computing  the  eigenvalue  X^  and  corresponding  eigenfunctions  e.  of  the  matrix 


where 


R = E {XXT} 


K - col  (x,  y,  z) 


(V-l) 


(V-2) 


is  a vector  of  east,  north,  and  vertical  displacem°nt  components,  and  E is  the  average  over  a 

time  window.  The  time  window  slides  down  the  record,  and  new  eigenvalues  and  eigenvectors 

are  computed  for  each  window.  At  a given  time,  the  eigenvector  e^  associated  with  the  largest 

eigenvalue  X^  is  the  direction  of  the  most  highly  polarized  wave,  which  can  be  P or  S motion. 

In  statistics,  e^  is  called  the  principal  component  and  has  variance  X ^ . If  X ^ » X 2 , then  the 

wave  is  highly  polarized.  A single-point,  adaptive  filter  which  passes  highly  polarized  waves 
4 

is  given  by  Madariaga  as 

r = i . :v-3> 

The  filtered  output,  components  are  calculated  by  first  projecting  the  principal  component 
onto  each  displacement  component,  and  then  scaling  each  component  by  f. 

Figure  V-3  shows  an  example  of  Iiis  program  as  applied  to  long-period  seismograms  of 
the  explosion  Carnikin,  6 November  1971,  recorded  at  subarray  1A  at  NORSAR.  Traces  1 to  3 
are  the  vertical,  easf,  and  . orth  components,  respectively,  of  the  P-wave.  After  rotations  in 
the  horizontal  and  vertical  planes,  the  radial  R,  vertical  Z and  transverse  T components  are 
displayed  in  traces  4 to  6.  These  components  are  orthogonal,  the  radial  component  defined  to 
be  toward  NORSAR  along  the  incident  ray  direction.  Traces  7 to  9 are  the  projections  of  the 
principal  component  e^  onto  the  rotated  data  using  a correlation  window  of  20  sec.  One  sees 
that  Cj  is  essentially  parallel  to  the  R component  fir  the  whole  length  of  the  record  shown, 
since  the  projections  onto  Z and  T components  are  nearly  zero.  Trace  10  is  the  fiRer  value 
f as  a function  of  time,  and  trace  11  is  the  product  of  f with  the  radial  component  of  trace  7. 

These  results  show  that  the  entire  coda  is  P-wave  motion  polarized  in  the  radial  direction. 
The  program  was  used  to  determ. ne  the  dip  angle  of  the  long-period  first  motion,  assuming  a 
back  azimuth  of  7.38°  toward  the  epicenter  of  Cannikin,  based  on  short-period  arrivals.  All 
subarrav  <jites  except  2B  were  used,  yielding  a dip  of  73.8°  ± 6.5°  for  21  sites.  Three  f tes 
included  showed  anomalously  large,  distorted  vertical  components,  due  apparently  to  a non- 
linear response  of  the  instruments  to  strong  motion.  Removing  these  sites  gave  a dip  angle  of 
72.9  ± 3.5°.  The  small  standard  deviation  suggests  that  these  data  *rom  Cannikin  are  usable 
for  long-period  source  studies.  ^ w.  Frasier 


C.  SURFACE  WAVES  FROM  CENTRAL  ASIAN  EARTHQUAKES 

8 

Following  the  method  of  Weidner  and  Aki,  w-e  have  studied  long-period  Rayleigh  waves 
from  pairs  of  earthquakes  in  Central  Asia.  Our  goal  is  to  determine  the  source  parameters 
(depth,  slip,  dip,  strike)  of  these  earthquakes,  and  then  to  evaluate  the  method  by  comparing 
our  results  against  independent  observations. 
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The  method  makes  use  of  amplitude  and  phase  spectra  of  surface  waves  from  earthquakes 
which  are  located  dost-  together  and  which  are  known  to  have  different  fault  plane  solutions.  In 

Q 

particular,  the  differential  phase  at  a given  observation  point  is  defined  by  Weidner  and  Aki  as 
follows: 


A?  = - «...  + »(t,  - T2  - 

where 

<p i*  <^2  = observed  phase  delay  for  events  1 and  2 in  cycles/recond 

Tl'T?  “ difference  between  origin  time  and  start  of  digitization  for 
events  1 and  2 in  seconds 

xl,x2  ~ ePi°entral  distances  for  events  1 and  2 in  kilometers 

u - frequency  in  hertz 
c = phase  velocity  in  kilometers /second. 


By  taking  differential  phase  and  amplitude  ratios,  we  hope  to  eliminate  the  propagation  effects 

and  retain  the  charactemstics  of  the  two  sources.  Using  theoretical  calculations  based  on  the 
o 

work  of  IV1.  Saito,  we  search  parameter  space  for  a best  fit  in  the  least-squares  sense  lo  all 
the  observed  spectral  ratios. 


TABLE  V-l 

SOURCE  PARAMETERS  FOR  TWO  CENTRAL 

ASIAN  EVENTS 

Event 

Number 

Date 

Origit. 

Time 

Latitude 

Longitude 

mb 

21 

11  May  1967 

14:50:57 

39.3‘N 

73.7° E 

5.5 

22 

14  Sept  1969 

16:15:26 

39.7  ° N 

74.8°  E 

5.5 

We  report  the  progress  on  one  earthquake  pair  located  in  the  Pamir  and  south  Tien  Shan 
mountains.  Information  about  this  pair  of  events  is  given  in  Table  V-l.  The  fault  plane  solu- 
tions were  determined  by  Molnar,  et  aL*°  We  have  digitized  vertical  component  Rayleigh 
waves  from  these  events  at  22  WWSSN  stations,  the  majority  of  which  are  located  on  the  Euro- 
pean continent.  Since  the  records  suffered  from  severe  interference,  a time  variable  filter 
(landisman,  et  al.*  * ) was  designed  and  applied  to  each  seismogram  so  as  to  pass  only 
fundamental  mode  Rayleigh  waves  between  a 30-  to  100-sec  period.  Sp*  _tral  ratios  were  com- 
puted for  each  station  in  this  period  range,  and  the  theoretical  calculations  led  to  the  fit  shown 
in  Fig.  V-4.  Each  graph  is  produced  by  changing  a parameter,  whilt  holding  the  others  fixed 
at  their  best  value.  Arrows  on  this  graph  indicate  minimum  in  the  standard  deviation  for  dif- 
ferential phase.  We  may  note  here  that  USCGS  reports  a depth  of  21  km  for  event  21,  while 
short -period  LASA  records  for  event  22  reveal  a phase,  possibly  pP,  which  places  this  event 
roughly  30  km  deep. 

We  are  encouraged  by  our  results  to  date,  and  continue  with  closer  analysis  on  more  earth- 
quake  pairs.  II.  J.  Patton 
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(a) 


Fig.  V-2(a-c ).  Far-field  source  pulses  and  synthetic  seismograms 
based  on  theoretical  fault  model. 
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ORIGINAL  SEISMOGRAMS 


ROTATED  SflSMOGRAMS 


PRINCIPAL  COMPONENT 


r I LT E R * RADIAL  PR.NCIPAL  COMPONEKT 


Fig.  V-3.  Example  ( F polarization  method.  Traces  1 to  3 are  vertical,  east,  and  nor*... 
components  of  long-period  P-wave  from  explosion  Car.  likin,  6 November  1971,  recorded 
at  subarray  1A  of  NORSAR.  Rotated  components  are  traces  4 to  6,  R being  radial  com- 
ponent toward  receiver  along  ray  path;  T i3  transverse  horizontal  component,  and  Z is 
upward  in  vertical  plane,  orthogonal  to  T and  R.  Traces  to  9 are  projections  of  prin- 
cipal component  onto  R,  Z,  and  T axes.  Trace  10  is  T,ola  ization  rilter,  and  tra~e  11  is 
product  of  traces  7 and  10. 
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Fig.  V-4.  Results  of  fit  to  obser  ved  spectral  ratios.  o\  2 is  weighted  sum 
of  squares  over  frequency  and  station  of  bserved  log  amplitude  ratio  minus 
theoretical  log  amplitude  ratio.  is  s milar  weighted  sum  of  squares  of 

observed  differential  phases  minus  theoretical  differential  phase.  For  de- 
tails, these  quantities  are  given  in  We’*1*  ;r  and  Aki.8 
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VI.  GLOBAL  SE’SMIC  ACTIVITY 


A.  CHARACTERISTICS  OF  SHORT-TERM  VARIATIONS  IN  SE1S1V  IC  ACTIVITY 

It  is  well  known  that  seismic  activity  varies  from  time  to  time,  although  the  characteristics 
of  this  variation  are  not  well  understood.  The  presence  of  this  variation  impacts  on  the  discrim- 
ination problem  in  several  ways.  In  order  to  evaluate  the  effectiveness  of  a seismic  network,  it 
is  common  to  count  the*  events  detected  during  a.  limited  period  (for  example,  the  International 
Seismic  Month1).  Estimates  of  the  level  of  background  seismic  activity  and  the  rate  of  occur- 
rence of  anomalous  events  are  frequently  obtained  in  the  same  way.  It  is  usually  not  clear 
whether  the  sample  discussed  is  a good  representation  of  the  seismic  activity  to  be  expected  at 
some  future  time.  It  is  possible,  for  example,  that  a network  processing  scheme  designed  to 
handle  the  long-term  average  rate  of  seismicity  in  real  time  may  become  swamped  with  events 
during  an  activity  maximum. 

Evidence2  has  been  given  that  the  PDE  earthquake  listing  is  approximately  complete  i'or 
magnitudes  (m^)  greater  than  5.1.  We  have  therefore  used  this  list  as  a basic  data  set.  Fig- 
ure Vl-1  shows  the  numb<  r of  events  listed  during  successive  40-dr.y  intervals  for  the  perioi 
1966  to  1974,  for  m^  > 5.2,  m^  > 5.5,  and  m,}  > 5.T.  The  period  of  40  days  was  chosen  as  being 
short  enough  to  show  rapid  fluctuations  in  activity,  but  long  enough  to  include  sufficient  events  for 
an  adequate  data  sample. 

Figure  Vl-1  has  many  features  ot  interest.  First,  the  fluctuations  in  event  numbers  are 
conside  pile.  For  example,  events  with  m^  > 5.5  range  from  a minimum  of  9 to  a n.aximum 
of  66  nur  40-oay  period.  Second,  tne  changes  are  frequently  very  abrupt.  Particularly  rapid 
changes  are  shown  during  Apr1”  1968,  June  1971,  July  1972,  and  June  1?7J.  Potential  reasons 
for  these  changes  are  under  investigation. 

Perhaps  the  most  important  aspect  of  Fig.  VI-1  is  that  the  fluctuations  in  activity  appear 

approximately  independent  of  magnitude.  To  illustrate  this.  Fig.  VI-2  shows  the  ratio  of  the 

number  of  events  with  m,  > 5.2  to  the  number  with  m,  > 5.5,  and  also  the  ratio  for  m,  > 5.5  to 
b b b 

m^  >5.8  With  the  exception  of  a few  unusually  large  values,  these  actios  s.h»  w renr  rkable 
stability.  This  implies  -hat  the  frequency-magnitude  curve  is  fluctuating  rap  dly  parallel  to  the 
frequency  axis,  with  very  little  change  in  shape.  Vi  he  re  is  some  evidence  for  a long-term 
change  in  shape,  which  is  small  but  persistent.)  We  prefer  to  use  the  term  "change  in  shape" 
rather  than  change  in  b-value  s\nce  we  have  shown  e rFer2  that  the  global  frequency-magnitude 
(m^)  curve  is  not  linear. 

We  conclude  from  this  study  tnat  seismic  activity  appears  to  have  real  sh.^rt-term  varia- 
tions in  level  that  are  considerable  in  size.  Further  investigations  of  this  phenomenon  are 


continuing. 


M.  A.  Chinnery 


P.  SPECTRAL  ANALYSIS  OF  EARTHQUAKE  OCCURRENCE  RATES 

In  the  previous  SATS,2  the  level  of  seismic  acuity  was  shown  to  be  a time-deperdeni  func- 
tion, independent  of  the  magnitude  range  over  which  it  is  determined.  Evidence  showing  that 
periodicities  exist  in  this  function  would  be  of  fundamental  importance  to  several  seismic  dis- 
crimination problems.  For  example,  the  ability  to  hide  an  explosion  in  an  earthquake  is  highly 
dependent  on  the  rate  at  which  suitably  large  enough  earthquakes  occur.  Knowledgi  of  a time 
interval  during  wnich  large  earthquakes  are  highly  likely  t<  occur  would  make  such  a scheme 
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more  feasible.  Evidence  for  periodicities  in  the  rate  function  would  also  indicate  what  controls 
the  initiation  of  rupture.  If  the  rate  function  is  random,  one  can  conjecture  that  the  physical 
characteristics  in  the  focal  region  control  whether  or  not  an  earthquake  will  occur.  On  the  other 
hand,  if  it  is  periodic  then  the  occurrence  of  earthquakes  may  be  more  directly  in  response  to 
regional  stress  fields  generated  perhaps  by  flow  in  the  asther.osphere.  Of  course,  the  real  sit* 
uation  is  at  best  a mixture,  the  proportions  of  the  mix  being  determined  by  the  relative  power 
contained  in  the  periodic  part  of  the  spect’um  of  the  rate  function. 

The  problem  of  determining  whether  or  not  a function  is  random  or  periodic  can  be  ap- 
proached by  comparing  an  estimate  of  power  at  a given  period  with  the  power  that  might  by 
chance  appear  at  that  period  if  an  equal  length  of  random  data  was  analyzed.  The  spectrum  from 
a finite  segment  of  random  data  will  be  peaked  and  not  necessarily  flat.  The  fact  that  it  is  ran- 
dom just  implies  that  the  average  of  a large  ensemble  of  spectra  of  the  sequence  will  tend  to  be 
white.  Consequently,  w.ien  dealing  with  finite  lengths  of  data,  spectral  peaks  must  be  assitmed 
a degree  of  significance  compared  with  a specific  random  process.  Many  authors  have  pointed 
out  that  spectral  peaks  for  earthquake  occurrence  rate  data  are  highly  significant  compared  with 
those  expected  from  random  data,  with  a Poisson  distribution  function  assuming  independent 
events.  They  usually  conclude,  though,  that  if  the  events  were  somev/hat  mutually  dependent, 
the  periodicities  would  not  be  significantly  nonrandom.  The  difficulty  in  determining  that  random 
process  to  compare  data  can  be  circumvented  somewhat  if  several  estimates  of  the  spectrum  of 
the  data  can  be  made  from  partially  or  totally  independent  se  gments.  The  nonrandom  spectral 
peaks  should  be  present  in  each  estimate.  Peaks  from  random  data  should  not  coincide,  and 
simultaneously  have  the  same  shape  so  that  the  average  spectrum  will  tend  to  be  white.  This 
behavior  Is  independent  of  the  probability  distribution  function  of  the  random  process. 

The  practical  problem  is  that  if  the  length  of  the  available  rate  data  is  short,  then  a reli- 
able spectral  estimate  is  difficult,  if  not  impossible,  to  make.  A recent  innovation,  due  to 
John  Burg  of  Stanford  University,  for  coriputing  the  autocorrelation  func  ion  for  short  sequences 
makes  it  possible  to  overcome  this  problem.  The  following  figures  present  the  results  of  using 
Maximum  Entropy  Spectral  Analysis,  using  the  Burg  technique  for  computing  the  prediction  erro»~ 
wavelet,  to  determine  the  nonrandom  periodic  behavior  of  Global  Earthquake  Occurrence  Kate 
data. 

Figure  VI-3  shows  the  earthquake  occurrence  rate  computed  using  all  earthquakes  in  the 
PDE  with  m^  ^ 5.5  from  1 ">64  to  the  present,  m^  ^ 5.5  was  chosen  because  the  catalog  is  es- 
sentially complete  above  this  threshold.  Taking  a limited  range  is  justified  by  evidence  that  the 
time  dependence  of  the  frequency  of  earthquakes  is  independent  of  magnitude  range.  We  chose 
to  start  at  1966  since  a large  single  swarm  in  the  Kat  Islands  dominates  the  data  previous  to 
this  date.  The  rate  was  determined  usirg  a time  interval  of  40  days.  The  total  number  of  points 
amounts  to  72  ard  (rovers  about  8 years.  Forty-day  averagin  * "»lso  acts  a.  an  anti-aliasing  filter 
for  any  high-frequency  power  that  might  exist  at  the  tidal  perieds.  The  idea  was  to  divide  these 
72  points  into  several  segments,  make  a spectral  estimate  of  each,  and  look  for  the  coincidence 
of  spectral  peaks.  If  the  rate  function  is  random,  the  spectral  average  will  tend  to  be  white. 
Figure  VI-4  shows  a random  sequence  of  numbers  with  a Poisson  probability  distribution  func- 
tion that  was  treated  exactly  the  same  as  the  rate  data  to  illustrate  the  previously  mentioned 
behavior. 

Figure  VI-5  illustrates  why  one  needs  to  make  ^veral  estimates  of  the  spectra  and  why 
simple  FFT  estimates  a.*e  lot  sufficient  to  do  th.s.  The  upoer  curve  is  the  FFT  spectrum  based 
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on  64  points  of  the  rate  data.  Directly  below  it  is  the  Burg  Maximum  Kntropy  (HMD)  spectrum 
for  the  same  data.  At  the  bottom  is  the  HMD  spectrum  for  64  points  of  the  random  data.  Ky 
comparing  the  bottom  two  spectra,  it  is  clear  that  one  cannot  say  whether  the  apparent  perio- 
dicities in  the  rate  dafa  are  random  or  real.  The  FFT  shows  essentially  the  same  peaks  as  the 
Burg  spectrum,  but  with  much  less  resolution.  If  the  data  were  divided  into  two  32-point  seg- 
ments and  the  spectra  .vere  computed,  the  FFT  spectra  would  have  little  resolution  and  further, 
and  probably  worse,  the  end  c mditions  implicit  in  the  FFT  would  make  the  resulting  spectra 
unreliable.  There  is,  for  ex  miple,  only  one  point  in  the  FFT  spectrum  between  the  peaks  near 

0.1  and  0.13  cycles  per*  month.  In  the  HMD  spectra,  about  100  points  are  plotted  between  the 
two  peaks. 

Figure  YT-6(a-e)  shows  how  the  spectra,  computed  from  3 segments  of  the  random  data, 
do  not  coincide.  Fach  segment  is  40  points  long  so  that  the  middle  segment  overlaps  half  of  the 
first  and  last  segments.  A filte.  ength  of  25  was  used.  The  spectr  a are  clearly  uncorrelated. 
Fsing  exactly  the  same  procedure  for  dividing  the  rate  data  and  using  the  same  filter*  length, 
we  computed  the  spectra  shown  in  Fig.  YT-7(a-e).  There  is  a triple  coincidence  of  peaKS  foe 
periodicities  at  0.3  and  7.5  months.  The  other  peaks  have  no  nonrandom  significance.  The 
slight  shift  in  the  bottom  to  7.6  months  is  not  significant,  In  the  HMD  method,  when  noise  is 
present  and  relatively  long  filter's  are  used,  such  shifts  are  common,  but  not  serious  in  this 
study. 

Since  resolution  is  so  sharp  and  because  of  noise  and  filter  length  effects,  absolute  atnpli- 
tu  »es  using  the  HMD  method  are  often  unreliable.  However,  the  integrated  power*  under  each 
peak  proves  to  be  quite  accurate.  Figure  VI  8(a-b)  compares  the  accumulated  power*  of  the 
spectra  for  the  rate  data  and  the  random  data,  and  shows  that  about  30  percent  of  the  total  power 
is  contained  in  the  bands  9. 3 ± 0.3,  and  7.5  ± 0.3  months.  The  average  of  the  accumulated  power 
for  the  random  data  tends  to  a straight  line,  indicating  whiteness. 

In  summary,  to  a high  <lesree  of  certainty,  periodn  ities  at  7.5  and  T 3 mouths  exist  in  the 
earthquake  rate  from  1066  to  the  present,  and  these  periods  contain  approximately  30  percent 
of  the  total  power  in  the  data.  The  physical  significance  of  this  is  yet  to  be  determined,  but 
there  are  no  externally  produced  strains  in  the  crust  with  these  periods  and  it  seems  likely  that 
we  are  ol  serving  a characteristic  of  a global  internally  induced  strain. 

T.  F.  landers 
M.  A.  C’hinnerv 
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Fig.  Vl-1.  Farthquake  counts  for  successive  40-day  intervals  fiom  PDF 
earthquake  listing  for  period  January  196A  to  April  1974. 
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Fig.  VI-2.  Ratios  of  earthquake  counts  shown  in  Fig.  VI  i. 


Fig.  VI- 3.  KarthquaKf*  occur lence  rate  per  40-day  interval, 
1964  to  present. 
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Fig.  VI-4.  Random  numbers  with  mean  of  30  and  Poisson  probaidHty 
distribution  at  40-dny  internals,  1964  to  present. 


Fig.  VI-5,  (a)  F FT  power  spectrum 

of  54  points  of  earthquake  occurrence 
rate;  (b)  BMK  power  spectrum  of  same 
data;  and  (c)  BMP)  power  spectrum  of 
64  points  of  Poisson  random  data. 
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Fig..  VI-6.  BME  power  spectra  for 
three  40-point  segments  of  random 
data,  (a)  First  half,  (b)  middle,  and 
(c)  1? st  half  of  data. 


Fig.  VI-7.  BME  power  jpectra  for 
three  40-point  segments  of  earth- 
quake occurrence  rate.  Peaks  at 
9.3  and  .5  months  coincide,  indi- 
cating continuing  power  at  those 
periods  over  complete  time  span, 
(a)  First  half,  (b)  middle,  and 
(c)  last  half  of  data. 
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Fig.  VI-8.  Integrated  power  from  0 to  f of  BME  spectra  for  occurrence  i ate 
and  random  data.  Approximately  30  percent  of  total  power  in  occurrence  rate 
lies  in  bands  9. 3 ± 0.3,  and  7.5  ± 0 3 months.  (a)  Earthquake  catalog  and 
(b)  Poisson  random  data. 


